Creation of gap solitons in Bose-Einstein condensates 
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We discuss a method to launch gap soliton-like structures 
in atomic Bose-Einstein condensates confined in optical traps. 
Bright vector solitons consisting of a superposition of two hy- 
perfine Zeeman sublevels can be created for both attractive 
and repulsive interactions between the atoms. Their forma- 
tion relies on the dynamics of the atomic internal ground 
states in two far-off resonant counterpropagating a + — de- 
polarized laser beams which form the optical trap. Numerical 
simulations show that these solitons can be prepared from a 
one-component state provided with an initial velocity. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj 



I. INTRODUCTION 

The Gross-Pitaevskii equation (GPE) has been used 
successfully in the recent past to explain various experi- 
ments on atomic Bose-Einstein condensates (see, e.g., the 
references in Jl|,||), and its validity for the description of 
the condensate dynamics at zero temperature is now well 
accepted. A further confirmation would be provided by 
the observation of solitary matter waves, the existence 
of which is generic to nonlinear Schrodinger wave equa- 
tions such as the GPE || . Such solitary waves could also 
find applications in the future, e.g., in the diffractionless 
transport of condensates. 

Various theoretical studies of this problem have al- 
ready been performed, predicting in particular the ex- 
istence of bright solitons, with corresponding spatially 
localized atomic density profiles, for condensates with at- 
tractive interactions B. Research on condensates with 
repulsive interactions has focused on the formation of 
gray solitons which correspond to dips in the atomic 
density. Their creation was investigated in Refs. 
their general properties were discussed in H, and Ref. 
Q worked out their analogy to the Josephson effect. 

Complementary and previous to this work, the forma- 
tion of atomic solitons was also examined theoretically 
in the context of nonlinear atom optics P-pT]|. In these 
studies the interaction between the atoms was assumed 
to result from laser-induced dipole-dipole forces, but this 
theory has not been experimentally tested so far. 

The reliance on attractive interactions to achieve 
bright matter-wave solitons in Bose-Einstein condensates 
is of course a serious limitation, due to the difficulties as- 
sociated with achieving condensation in the first place for 
such interactions. The purpose of the present article is 
the theoretical exposition of an experimentally realizable 
geometry that allows one to create bright gap soliton-like 



structures in Bose condensates, for both attractive and 
repulsive signs of the two-body scattering length. Gap 
solitons result from the balance of nonlinearity and the 
effective linear dispersion of a coupled system, e.g., coun- 
terpropagating waves in a grating structure, and appear 
in the gaps associated with avoided crossings. Gap soli- 
tons have previously been studied in a variety of phys- 
ical contexts, but particularly in nonlinear optics [ p2| . 
They were also studied in the framework of nonlinear 
atom optics 0, but in this case the two states involved 
are connected by an optical transition, and the effects of 
spontaneous emission can cause significant problems j^]. 

Several main reasons motivate our renewed interest in 
this problem. First, we already mentioned that bright 
gap solitons are known to exist in nonlinear systems irre- 
spective of whether the nonlinear interaction is repulsive 
or attractive . With regard to atomic condensates this 
means that they should be observable, at least in prin- 
ciple, also for Na and Rb where the positive interatomic 
scattering length gives rise to a repulsive mean interac- 
tion. Further, the study of bright solitary waves is of 
interest as they might be easier to detect than gray ones, 
and they could find future applications, e.g. in atomic in- 
terferometry [0. An additional reason to study atomic 
gap solitons is the fact that they consist inherently of 
a superposition of two internal states, in our case two 
different Zeeman sublevels of the atomic ground state. 
As such, they offer a further example of a multicompo- 
nent Bose condensate the study of which has already re- 
ceived much interest recently [[14-16 . Finally, the recent 
demonstration of far-off resonant dipole traps for con- 
densates opens up the way to the "easy" generation and 
manipulation of such spinor systems. 

This paper is organized as follows. Section II describes 
our model. The physics relevant for the generation of gap 
solitons as well as orders of magnitudes for the various 
experimental parameters involved are discussed in Sec. 
Ill, while Sec. IV presents a summary of our numerical 
results. Finally, conclusions are given in Sec. V. 



II. THE MODEL 

The situation we consider for the generation of atomic 
gap solitons makes use of the recently achieved con- 
finement of Bose condensates in far off-resonant optical 
dipole traps |l7]]. We consider explicitly a trap consist- 
ing of two focused laser beams of frequency w; counter- 
propagating in the Z-direction and with polarizations o~ + 
and <j~ , respectively. These lasers are used to confine a 
Bose condensate which is assumed for concreteness to 
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consist of Na atoms. The condensate is initially pre- 



pared in the \g, F g = l,M g 



1) atomic ground state. 



For lasers far detuned from the resonance frequency cu a 
of the nearest transition to an excited hyperfine multiplct 
|e, F e , M e — —F e , . . . , F e ) the dynamics of a single atom 
in the trap can be described by an effective Hamiltonian 
of the form El 



H. 
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which is derived by adiabatically eliminating the excited 
states in the dipole and rotating wave approximations. 
In the Hamiltonian (|l|), the operators R and P denote 
the center-of-mass position and momentum of the atom 
of mass m, the ket \j) labels the magnetic sublevel of 



the Na ground states, \j) 
Ki = uji/c. Furthermore, 
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<5'(R) = «5s(R)/2, 



and 
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where we have introduced the detuning 5 = u>i — to a and 
the position-dependent saturation parameter 



S (R) = 



V 2 £ 2 (R) V 2 £ 2 (K) 



5 2 + T 2 /4 



S 2 



(3) 



In this expression, T> denotes the reduced dipole moment 
between the states \g) and |e), T is the upper to lower 
state spontaneous emission rate, and £(R) is the slowly 
varying laser field amplitude at point R, the plane-wave 
factors exp[i(±KiZ — uit)] having already been removed 
from the counterpropagating waves. In the following, we 
assume that £, which is identical for both fields, varies 
only in the transverse X- and ^-directions and is con- 
stant along the trap axis Z: This approximation is valid 
if the longitudinal extension of the confined BEC is much 
less than the Rayleigh range of the trapping fields, a con- 
dition we assume is satisfied. The numerical coefficients 
dj, which depend on the specific value of F e , are of the 
order of or somewhat less than unity. Note that except 
insofar as T appears in the saturation parameter s(R), 
the effects of spontaneous emission are neglected in this 
description^]. 



x In the discussion of the atomic dynamics and Eq. (|l|) we 
have assumed that the initial state is coupled only to one ex- 
cited hyperfine multiplet. However, in the optical trap the 
detuning of the laser frequency is large compared even to the 
fine structure splitting of the excited states, so that in princi- 
ple several different hyperfine multiplets should be taken into 
account. Fortunately, the coupling to any of these multiplets 
gives rise to addititional contributions to Eq. (^) which are of 
the same analytical structure as the one given above. Only 



The first term in the single-particle Hamiltonian ([[]) 
describes the quantized center-of-mass atomic motion, 
the second and third terms the (position dependent) 
light-shifts of the j = 0,±1 states, and the final term, 
proportional to d 2 , describes the coupling between the 
j = ±1 states by the counterpropagating fields. For ex- 
ample, coupling between the j — — 1 and j ' — 1 states 
arises from the process involving absorption of one cr + 
photon and subsequent re-emission of a a~ photon. How- 
ever, since the circularly polarized fields are counterprop- 
agating this process also involves a transfer of linear mo- 
mentum 2hKi along the Z-axis, and this accounts for 
the appearance of the spatially periodic factors, or grat- 
ings, exp(±2iKiZ) in the coupling terms. As shown be- 
low, these gratings provide the effective linear dispersion 
which allows for gap solitons in combination with the 
nonlinearity due to many-body effects. 

To describe the dynamics of the Bose condensate 
we introduce the macroscopic wave function ^(R, t) = 
[$i(R, t), fy-i(R., t)] T normalized to the total number of 
particles N, Here ty is omitted as it is coupled to ^/±i 
neither by H e ff n or by the nonlinearity if it vanishes 
initially, which we assume in the following. The time 
evolution of the spinor ^(R, t) is determined by the two- 
component Gross-Pitaevskii equation 

a* 

ffe//*(R,i) (4) 
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[U a \fh. (R, t)\ 2 + U b \^ 1 (R,t)\ 2 }i^ 1 (R, t) 
[Wi(R, t)\ 2 + <7 a |V-i(R, t)| 2 ]^-i(R, t) 



In the following we approximate the nonlinearity coeffi- 
cients by U a ~ Ub = U = 4:Trh 2 a S c/fn with a sc the s-wave 
scattering length. 

To identify the key physical parameters for gap soliton 
formation, and to facilitate numerical simulations, it is 
convenient to re-express Eq. (^) in a dimensionless form 
by introducing scaled variables r = t/t c , r = R/Z c and 
ipj — VP j/ ' yfpc with 

t c = l/(d 2 5' ), (5) 
l c = t c ■ KK/m, (6) 
Pc = \d 2 hS' /U\,, (7) 

where S' — <5'(R = 0). Note that for our choice of cfe = 
— 1/4, and for red detuning, we have > 0. Equation 
(0) then reads 
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-MA+MM e 2 ^5>(r)/S> 
- 2ik > z 5'(r)/6' -MA + ^P- 



the values of dj, T>, and F are different. This means that Eq. 
(1) may still be used in this case, the effects of the additional 
multiplets being included as modifications of the values of the 
coefficients dj. For simplicity, however, we will use the val- 



d = 1/2, di = 1/4, and d 2 



T/4 @. 
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+ [ S gn(d 2 K/[/)(|^i| 2 + |^i| 2 )] ) , (8) 

where A is the Laplacian in scaled variables, and we have 
introduced the dimcnsionless mass-related parameter 



M = d 2 S' Q m/(2hKf), 
so that h = K t l c = 1/(2M). 



(9) 



III. GAP SOLITONS 



exact gap solitons solutions are coherent superpositions 
of the j = ±1 states where the phase and amplitude of 
the superposition varies spatially in a specific manner: it 
is not a priori clear that these gap solitons can be ac- 
cessed from an initial state purely in the j = — 1 state for 
example. Furthermore, inclusion of transverse variations 
and spatial derivatives beyond the slowly-varying enve- 
lope approximation introduced above could, in principle, 
destroy the solitons |19|. For the numerical simulations 
to be presented here we work directly with Eqs. (^) which 
does not invoke these approximations. 



In this section we discuss the conditions under which 
Eqs. (||) yield gap soliton solutions. Rather than re- 
producing the explicit analytic forms of these solutions, 
which are readily available in the literature, here we in- 
troduce the reduced equations which yield gap solitons, 
and discuss the physics underlying their formation. Esti- 
mates for the orders of magnitude of various parameters 
characterizing atomic gap solitons are also given. 



A. Reduced soliton equations 

Two key approximations underly the appearance of 
gap solitons: First, we neglect all transverse variations 
of the electromagnetic and atomic fields, thereby reduc- 
ing the problem to one spatial variable z. Furthermore, 
we can set S'(r)/S' = 1. Second we express the atomic 
fields in the form 

4>±i(z,t) = exp{i[±fci* - (1/(4M) - 1)t]}0 ±1 (z,£), 

(10) 

and we assume that the atomic field envelopes tp±i(z, t) 
vary slowly in space in comparison to the plane- wave fac- 
tors that have been separated out, so that only first-order 
spatial derivatives of the field envelopes need be retained 
and only the spatial harmonics indicated included. Un- 
der these assumptions Eqs. (@) reduce to 
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+sgn(d 2 ?^/[/)(|0i| 2 + | ( M 2 ) 



(11) 



Aceves and Wabnitz 13 have shown that these dimen- 
sionless equations have explicit travelling solitary wave 
solutions of hyperbolic secant form. Thus, the optical 
trapping geometry we propose here can support atomic 
gap solitons under the appropriate conditions. 

Having established that our system can support gap 
solitons our goal in the remainder of this paper is to 
demonstrate through numerical simulations that these 
solitons, or at least a remnant of them, can arise for real- 
istic atomic properties and that they can be created from 
physically reasonable initial conditions. In particular, the 



B. Intuitive soliton picture 

A simple and intuitively appealing explanation of the 
reason why Eq. (jg) supports soliton solutions goes as 
follows: Consider first the one-dimensional nonlinear 
Schrodinger equation 



-Md 2 ip/dz 2 + g\ip\ 2 ^J. 



(12) 



This equation has bright soliton solutions if the effects of 
dispersion and nonlinearity can cancel each other. For 
this to happen, it is necessary that Mg < 0. In the usual 
case the mass-related coefficient M is positive, so that 
bright solitons can only exist in condensates with attrac- 
tive interactions g < 0. But consider now the dispersion 
relation u(k) = Mk 2 ± vl + k 2 for the linear part of Eq. 
(||) obtained after neglecting the transverse dimensions 
and performing the transformation 



ip±i — a±i exp{i[fc±iz — uj(k)T}}. 



(13) 



Thereby, k±\ = k ± ki for the j = ±1 states, and A; is a 
relative longitudinal wave vector. The dispersion relation 
consists of two branches, which in the absence of linear 
coupling take the form of two parabolas corresponding 
to the free dynamics of the internal states | ± 1). How- 
ever, the linear coupling between these states results in 
an avoided crossing at k = 0, see Fig. |l|. If the system is 
in a superposition of eigenstates pertaining to the lower 
branch of the dispersion relation, then at the crossing it 
can be ascribed a negative effective mass. One can thus 
expect that in this case the system can support soliton 
solutions even though the interaction is repulsive. For an 
attractive interaction, in contrast, soliton creation should 
be possible in all regions of the spectrum with positive 
effective mass. 

From the dispersion relation picture, one can easily 
infer further properties of repulsive interaction solitons. 
First, they will only exist for weak enough dispersion, as 
the lower branch of the dispersion curve has a region with 
negative curvature only as long as the dimensionless mass 
M < 0.5. Also, the maximum possible velocity can be 
estimated to be of the order of [1 — (2A/) 2 / 3 ] 3 / 2 which is 
the group velocity at the points ±k = ±\/ (2M)~ 2 / 3 — 1 
of vanishing curvature in the dispersion relation. Finally, 
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for a soliton at rest the contributions of the internal states 
|1) and | — 1) will approximately be equal, but solitons 
traveling with increasing positive, resp. negative, veloc- 
ity will be increasingly dominated by the |1), resp. | — 1), 
contribution. 

This qualitative discussion is in agreement with the 
analytic results of Ref. Jll|. More precisely, the solu- 
tions of Ref. |b| are solitary waves. In the following, we 
will be concerned with the creation of long-lived localized 
wave packet structures which are brought about by the 
interplay between nonlinearity and dispersion described 
above. We will continue to refer to these structures as 
gap solitons for simplicity. 



sion. From this value, it is possible to obtain a first es- 
timate of the transverse confinement of the condensate 
required in our two-dimensional model: We assume that 
the transverse spatial dependence of the atomic density 
can be modeled as the normalized ground state density 
g(X) of the harmonic trap potential 5'(X) — mui 2 X 2 /2 
]5|. The soliton density can hence be roughly estimated 
as p(R) = e{Z)e{AZ - Z)g(X)N/AZ, where 6(Z)is 
the Heaviside step function, N a typical total number 
of atoms in the soliton and AZ its length. From this 
condition, and using the typical values for iV and AZ 
previoulsy discussed leads to a lower limit for lo x in the 
range between 100 and 1000 Hz. Altogether, these vari- 
ous estimates are well within experimental reach. 



C. Soliton estimates 

We now turn to a discussion of the typical orders of 
magnitude which characterize the soliton solutions of Eq. 
(Q) . We note from the outset that the analytical solutions 
of Ref. |l3[ , as well as our numerical simulations, indicate 
that these characteristic scales can be directly inferred 
from the scale variables in Eqs. (U)-(@) which bring the 
Gross-Pitaevskii equations into dimensionless form. For 
example, the spatial extension of the scaled wave function 
ip, as well as the total norm J dv{\tp\\ 2 + |V>-i| 2 ) are of 
order unity.0 

In order to obtain estimates for the characteristic 
length, time and density introduced in Eqs. we 
use the parameter values of the Na experiment of Ref. 
JlTj as a guidance. For Sodium, T = 61 MHz, the sat- 
uration intensity I s = 6.2 mW/cm 2 and the resonance 
wavelength A a = 589 nm. Choosing the trap wavelength 
to be far red-detuned from this value, with A; = 985 nm, 
and a maximum laser intensity / as 1 kW/cm 2 one ob- 
tains a characteristic scale t c for the time evolution of 
the condensate of the order of 50 /xs. The characteristic 
length is obtained by multiplication with the recoil ve- 
locity iv-ec =1.8 cm/s, which yields l c rs 1 — 2 /im. This 
yields the dimensionless mass M w 0.1. Finally, the or- 
der of magnitude of the characteristic density p c = 10 14 
cm -3 , which means that a soliton typically contains of 
the order of p c l\ w 100 — 1000 atoms. These estimates 
are confirmed by our numerical simulations, which show 
that the typical extension of a soliton is several l c in the 
^-direction, about one l c in the transverse direction and 
it contains about 1000 atoms. 

Our numerical simulations show that the maximum 
dimensionless atomic density in a soliton is always of 
the order of p m = 0.1, which appears to produce the 
nonlinearity necessary to balance the effects of disper- 



2 Note that the precise values of di and are only of rele- 
vance for the scaling between Eqs. (JiJ) and (^). They do not 
influence the essential physics of the system. 



IV. NUMERICAL RESULTS 

Having characterized the idealized gap soliton solu- 
tions of Eq. (Q) we now investigate whether they can be 
accessed from realistic initial conditions. To this end, we 
study numerically the following situation. A condensate 
of Nq atoms in the internal state | — 1) is initially prepared 
in a conventional optical dipole trap which provides only 
a tranverse confinement potential V c . This potential is 
assumed to be Gaussian, with a trapping frequency u> x 
at the bottom. Axially, the condensate is confined by a 
harmonic magnetic trap of frequency uj z . At t = the 
magnetic trap is turned off and the polarizations of the 
trapping light fields are switched to the er + — o~ config- 
uration, with V c being unchanged. 

The simulations are performed in two spatial dimen- 
sions only, x and z, as this can already be expected to 
capture the relevant physics without requiring excessive 
computational resources. We concentrate on the more in- 
teresting case of a condensate with repulsive interactions 
since that is the new case in which solitons are expected. 
In order to estimate atom numbers the transformation 
between Eqs. (||) and (0) is performed after replacing U 
by U / y/ (X 2 ) where the variance (X 2 ) is determined from 
the two-dimensional wave packet structure in question. 

The main purpose of the numerical simulations is to 
show that gap solitons can be formed out of condensate 
wave functions whose initial parameters lie within a rel- 
atively broad range. It is only necessary to choose u) x , 
lo z , and No such that the spatial extension of the ini- 
tial condensate and the atom number are comparable to 
typical soliton values. They need not take on precisely 
defined values and the initial wave function does not have 
to match closely the form of a soliton. 

However, the condensate will not couple effectively to 
a soliton if it is at rest initially. Such a situation cor- 
responds to the point with k = ki in Fig. |l| where the 
effective mass is still positive (fc; > ko). The key to an 
efficient generation of solitons is therefore to provide the 
condensate with an initial velocity V, n = HKi n /m close 
to the recoil velocity V r in the Z-direction. This may 
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be achieved, e.g., by suddenly displacing the center of 
the magnetic trap. The initial wave function can then 
be written as ip-i(x, z,0) — ip g (x, z, 0) exp(iki n z) with 
ip g (x,z,0) the ground state of the combined optical and 
magnetic trap and fcj„ = Ki n l c Q. It is thus placed 
in the vicinity of the avoided crossing. Experimentally, 
condensates have already been accelerated to velocities 
in this range by a similar method in connection with the 
excitation of dipole oscillations [^(J . 

Figure || shows an illustrative example for the forma- 
tion of a soliton out of the initial distribution. It depicts 
the evolution of the transverse averaged atomic density 

N(z;t)= f dx(\Mx,z;r)\ 2 + \tP-i(x,z;t)\ 2 ), (14) 

as a function of the scaled variables z and r. In this ex- 
ample, A; = 985 nm, the maximum intensity / = 0.88 
kW/cm 2 , uj x = 6000 s-\ lo z = 5100 s~\ and V m = 
0.9V r . The characteristic scales are thus l c = 1.4 /jm 
and t c = 77 /xs, the coefficient M = 0.06. We choose 
J dxdz\tp-i(x, z, 0)| 2 = 6.21 which corresponds to an 
initial atom number of 2900, approximately. Figure ^ 
shows the formation of a soliton after an initial tran- 
sient phase having a duration of 50 t c hari approximately. 
This transient phase is characterized by strong "radia- 
tion losses" . They occur because half of the initial state 
pertains to the upper branch of the dispersion relation 
u>(k) which cannot sustain solitons. The shape of the 
created soliton is not stationary in time but appears to 
oscillate. Further examination shows that norm of the 
j = — 1 state is slightly larger than that of the j = 1 
state, as is expected for a soliton moving slowly in neg- 
ative z-direction |D||. Figure || shows the atomic den- 
sity P(x, z) = \ipijx, z; t)\ 2 + \ip_i(x, z; r)| 2 in the soli- 
ton at r = 100. The soliton contains about 500 atoms. 
The inset depicts the longitudinally integrated density 
/ dz(\ipi(x, z; r)| 2 + \ip_i(x, z; r)| 2 ) and the transverse 
confinement potential in the shape of an inverted Gaus- 
sian. The soliton thus spreads out over half the width of 
the potential well, approximately. 

Various numerical simulations were performed in order 
to assess the dependence of soliton formation on the var- 
ious initial parameters. When changing the atom num- 
ber iVo a soliton was formed over the whole investigated 
range between 500 and 4000 atoms. With increasing No 
and thus increasing effects of the nonlinearity the final 
velocity of the soliton changed from negative to positive 
values. For large No a tendency to form two soliton wave 
packets out of the initial state was observed, however, 
the formation of each of these solitons is accompanied by 
large radiation losses which destabilize the other one. As 
to the transverse confinement parameter tu x soliton for- 
mation was observed for uj x > 4000 s _1 . At ui x — 2000 
s _1 a stable structure was no longer attained which is in 
rough agreement with the estimate given above. Whereas 
these numbers indicate a relatively large freedom in the 
choice of No and oj x (for lo z similar results can be ex- 
pected) a somewhat more restrictive condition is placed 



on the initial velocity Vi n . Its value should be chosen 
from the interval between 0.8 and 1.0 V r in order to 
guarantee soliton formation, the lower bound being de- 
termined by the point of vanishing curvature in the dis- 
persion relation. For Vi n > V r the initial wave function 
is situated more and more on the upper branch of the 
dispersion relation so that the tendency to form solitons 
is diminished rapidly. 



V. SUMMARY AND CONCLUSION 

In conclusion, we have demonstrated that gap soliton- 
like structures can be created in a Bose condensate con- 
fined in an optical dipole trap formed by two counterpop- 
agating a + — cr~-polarized laser beams. Bright solitons 
can be formed not only for atomic species with attractive 
interactions but also in the repulsive case. This is ren- 
dered possible because the atoms can be ascribed a neg- 
ative effective mass if their velocity is close to the recoil 
velocity. The repulsive interaction solitons are inherently 
superpositions of two hyperfine Zeeman sublevels. The 
discussion of characteristic scales and numerical simula- 
tions indicated that the actual observation of these struc- 
tures should be achievable within the realm of current 
experimental possibilities. 

In our theoretical treatment spontaneous emission was 
neglected, an approximation justified by the large detun- 
ings in the optical trap ^] . The effects of anti-resonant 
terms, which were also ignored, might be of more im- 
portance. This question, as well as three-dimensional 
numerical studies, could be the subject of future work. 
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FIG. 1. Dispersion relation of the linear part of Eq. Q for 
the parameter value M = 0.06. The fc-value associated with 
a condensate at rest is denoted by fc;, the point of vanishing 
curvature is indicated by fco- 

FIG. 2. Formation of a soliton out of an ini- 
tial distribution. Depicted is the integrated density 
N(z; t) = J dx(\ip-i_(x, z; r)| 2 + \ip-i(x, z; r)| 2 ) in scaled vari- 
ables. Parameter values given in the text. 



FIG. 3. Atomic 

density P(x^z) — \ip\(x, z; t)| 2 + \ip-i(x, z; r)| 2 of the soli- 
ton of Fig. H at r = 100. The inset shows the integrated 
density Q(x) = J dz(\ipi(x, z; r)| 2 + \ip-i (x, z; r)| 2 ) and the 
transverse confinement potential V(x). 
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